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We report the experimental implementation of Grover's 
quantum search algorithm on a quantum computer with three 
quantum bits. The computer consists of molecules of ^Re- 
labeled CHFBr2, in which the three weakly coupled spin-1/2 
nuclei behave as the bits and are initialized, manipulated, and 
read out using magnetic resonance techniques. This quan- 
tum computation is made possible by the introduction of two 
techniques which significantly reduce the complexity of the 
experiment and by the surprising degree of cancellation of 
systematic errors which have previously limited the total pos- 
sible number of quantum gates. 

Elementary quantum computations K^ have recently 
been demonstrated experimentally Q using coupled nu- 
clear spins as quantum bits (qubits) and solution nuclear 
magnetic resonance (NMR) techniques to prepare, ma- 
nipulate, and detect the spins |H,|j. Aside from the well- 
understood scaling limitations due to the use of a high- 
temperature (almost random) system instead of a low- 
temperature (low entropy) polarized spin system MM, 
the crucial limitation in applying this method to imple- 
ment larger quantum algorithms has been systematic er- 
rors in the quantum gates. These gates are implemented 
by applying pulses of radiofrequency (RF) electromag- 
netic fields of precise duration and phase, which are in 
practice highly inhomogeneous over the sample volume, 
causing the gate fidelity [g| to be less than 95%. Produc- 
ing a homogeneous field is difficult because of the sample 
geometry and the necessity of keeping the field transverse 
to the 11.8 Tesla alignment field. If such systematic er- 
rors simply accumulated, these observations would imply 
that fewer than 90 gates p applied to any one spin could 
ever be cascaded in these systems. 

One technique which has been proposed for control- 
ling errors in quantum gates is quantum error correction 
prt, and in fact, it is known that if the error proba- 
bility per gate is smaller than a certain threshold pi[ , 
and in addition the architecture of the computer satis- 
fies certain criteria, then quantum computation could be 
carried out robustly and indefinitely, despite the errors. 
Unfortunately, this threshold is currently estimated to be 
below 10""*, which is virtually unreachable with current 
NMR quantum computation techniques and present lev- 
els of systematic errors. Now, the principle attribute of 



traditional error correction techniques is their ability to 
correct completely random errors which originate from 
fundamentally irreversible decoherence phenomena [ [l2[ ; 
in principle, systematic errors, which are inherently re- 
versible — at least on an appropriate time scale — should 
be much easier to control, given knowledge about their 
origin. 

Here, we report on the experimental realization of a 
three-qubit quantum search algorithm, and the observa- 
tion of 28 full cycles of the theoretically predicted oscilla- 
tory behavior for the algorithm. This was made possible 
by the application of two techniques to reduce the num- 
ber of gates required for the implementation, and by the 
control and cancellation of systematic errors. More than 
280 two-qubit quantum gates involving 1350 RF pulses 
were successfully cascaded, which exceeds not only the 
number of gates used in all previous NMR quantum com- 
puting experiments but also the limitation of 90 gates, 
imposed by cumulative systematic errors. 

The experiments were carried out with molecules con- 
taining n — 3 weakly coupled spin-1/2 nuclei, subject to a 
strong magnetic field Bq- Each spin represents one qubit, 
with the ground and excited state serving as |0) and |1). 
The Hamiltonian is well approximated by (?i = 1) [p[ 
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where the first term represents the free precession of each 
spin i about —Bq, with I^i the angular momentum op- 
erator in the z direction and cui the Larmor frequency of 
spin i. The second term describes a scalar spin-spin cou- 
pling of Jij between spins i and j and TCenv represents 
coupling to the environment, which causes decoherence. 
A resonant RF field rotating in the x-y plane can be gated 
on to perform single-spin rotations. 

The first method to reduce the number of one- and two- 
spin gates used to realize any given n-qubit unitary oper- 
ation starts with a library of efficient implementations for 
often-used building blocks. For example, four equivalent 
realizations of a Z-rotation using an RF coil in the trans- 
verse plane, are Z = YXY = YXY = XYX = XYX, 
where X and X denote +90° and -90° rotations about the 
X axis. _And CNOT,y_= XjYjnjZ.Yj = XjYjTijZiYj = 
YjTijZiYjXj = YjTijZiYjXj where r^ describes a time 
evolution of 1/2 Jij, during which all couplings except 
Jij are refocused g (CNOTy represents a controUed- 
NOT which flips spin j if and only if the control spin i 
is |1)). For any given quantum circuit, e.g. Fig. |l| (af- 
ter jl^]), specific implementations of the building blocks 
can then be laid out such that adjacent pulses invert 
each other, so they can be omitted. The simplification 
process can exploit the commutation of Z rotations with 
time evolutions (the Hamiltonian) , and of RF pulses on 
different spins with each other. Furthermore, NMR per- 
mits a direct implementation of many useful gates such as 



the controUed-y (Fig. |]), which can be reahzed by Ty/2 
and a few single-spin rotations. If only CNOT's and 1- 
qubit gates were used to implement the controlled-y's 
(as in ||l^) and if no pulses were canceled out, 70.5 7r/2 
pulses and 8 evolutions of 1/2J would be required for 
the Toffoli gate (Fig. p. This was reduced to 19 pulses, 
2 evolutions of 1/2 J and 3 of 1/4J, by using the simpli- 
fication methodology 1 14 1 . 

The second method to reduce the complexity concerns 
the initialization of the qubits to the ground state, gen- 
erally the first step in quantum algorithms. The density 
matrix peg of the system of Eq. n^ in equilibrium at room 
temperature is highly mixed (hcui <ti ksT, with ksT the 
thermal energy). The use of room temperature nuclear 
spins is made possible by converting the equilibrium state 
into an "effective pure state" [^ (one population deviates 
from a uniform background). We propose the prepara- 
tion of an effective pure ground state pg by a new variant 
of temporal averaging [|l5| . The original scheme involves 
the summation of 2" — 1 density matrices, each of which 
is obtained from p^q by cyclicly permuting all popula- 
tions except the ground state population |l5[ ; but this 
gets rather involved for n > 3. A much simpler approach 
is to take a weighted sum of (diagonal) density matri- 
ces pi which are obtained from p^q by performing just a 
few CNOT's to rearrange the 2" — 1 populations. The 
weights wi are determined by solving a set of linear equa- 
tions ^i wi diag(p;) = diag(pg). This approach may still 
require up to 2" — 1 experiments to get exactly pg, but 
allows one also to approximate pg using far fewer exper- 
iments. All the data shown below were obtained using 
just three experiments, resulting in an expected variance 
of the 2" — 1 populations of only 7% of their average 
value. These modifications thus increase both the sim- 
plicity and scalability of temporal labeling. 

We used these methods to implement the eight in- 
stances of the 3-qubit Grover algorithm. This algo- 
rithm achieves a quadratic speed-up over classical algo- 
rithms when searching in an unsorted database with N 
entries |^: given an unknown function f{x) with the 
promise that /(xq) = — 1 for a unique xq, and f{x) — 1 
elsewhere, a classical search needs 0{N) (4.375 for N—8) 
attempts to find xq, while 0{^/N) (2 for N=8) queries 
suffice using Grover's quantum algorithm |3|. For A^=8, 
the algorithm requires four Toffoli gates |Q3|, plus several 
1-qubit gates. Earlier experiments [0M| required at most 
one Toffoli gate and a few two-qubit gates. 

We selected ^^C-labeled CHFBra |ll|| for our experi- 
ments. The Hamiltonian of the ^H-^^T^^^C spin system 
is of the form of Eq. |l| with Jhc = 224 Hz, Jhf = 50 Hz 
and Jpc — —311 Hz. The scalar interaction with the Br 
nuclei is averaged out and only contributes to Tienv Pi- 
Experiments were carried out at IBM using an 11.7 T 
Oxford Instruments magnet and a Varian Unity Inova 
spectrometer with a triple resonance (H-F-X) probe from 
Nalorac. The output state of each spin i, |0) or |1), can 



be determined from the phase of the signal induced in 
the RF coil after applying an Xi read-out pulse. In fact, 
the spectrum of the signal of any one spin suffices to de- 
termine the output state of all the spins given that they 
are in an effective pure energy eigenstate and that they 
are all mutually coupled: each spectrum then contains 
only a single line, the frequency of which, combined with 
the knowledge of the Jij , reveals the state of the remain- 
ing spins (Fig. || (a), inset). We also reconstructed the 
complete output deviation (traceless) density matrix M 
using quantum state tomography [^ (Fig. ^ (a)). 

The agreement between experimental results and the- 
oretical predictions is good, considering that about 100 
pulses were used and that the systematic error rate ex- 
ceeds 5% per RF pulse (the measured signal loss due to 
RF field inhomogeneity after applying Xi). This sug- 
gests that the systematic errors cancel each other out to 
some degree. We examined this in more detail in a series 
of experiments with increasingly longer pulse sequences 
executing up to 28 Grover iterations (repeated execu- 
tions of the two main steps of Grover's algorithm Q). 
Theoretically, the probability of \xo) oscillates as a func- 
tion of the number of iterations k, reaching a first max- 
imum for k = 0{^/N) §]. Fig. | (a) shows that the 
diagonal entry dx„ of pexp oscillates as predicted but the 
oscillation is damped as a result of errors, with a time 
constant T^ of 12.8 iterations. However, T^ would have 
been smaller than 1.5 if the errors due to just the RF 
field inhomogeneity were cumulative (Fig. (a), solid 
line). Remarkably, after a considerable initial loss, d^^ 
decays at a rate close to the ^^C T2 decay rate B (dashed 
line), which can be regarded as a lower bound on the 
overall error rate. A more complete measure to quan- 
tify the error and benchmark results is the relative error 
fir =11 cpexp-pth h / W Pthh, where pexp and pth are the 
experimental and theoretical (traceless) deviation density 
matrices, and || • II2 is the 2-norm ||l^. Comparison of 
er with c = 1 and c equal to the inverse of the signal loss 
(Fig. ^ (b)) reveals that signal loss dominates over other 
types of error. Furthermore, the small values of e^ with 
c > 1 suggest that \xq) can be unambiguously identified, 
even after almost 1350 pulses. This is confirmed by the 
density matrix measured after 28 iterations, which has a 
surprisingly good signature (Fig. (b)). Given the error 
of > 5% per single 7r/2 rotation, all these observations 
demonstrate that substantial cancellation of errors took 
place in our experiments. 

The error cancellation achieved was partly due to a 
judicious choice of the phases of the refocusing pulses, 
but a detailed mathematical description in terms of er- 
ror operators is needed to fully exploit this effect in ar- 
bitrary pulse sequences. This difficult undertaking is 
made worthwhile by our observations. This conclusion 
is strengthened by a similar observation in [Q. Also, we 
believe that error cancellation behavior is not a property 
of the Grover iterations, because we found experimen- 



tally that the choice of implementation of the building 
blocks dramatically affects the cancellation effectiveness. 
Whereas the cancellation of systematic errors makes it 
possible to perform surprisingly many operations, the 
methods for simplifying pulse sequences reduce the num- 
ber of operations needed to implement a given quantum 
circuit. This combination permitted the succesful realiza- 
tion of Grover's algorithm with 3 spins and brings many 
other interesting quantum computing experiments within 
reach. 
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FIG. 1. Quantum circuit for a three-qubit Toffoli gate 
which flips the phase of the 1 1) 1 1) j 1) term. Vertical lines con- 
necting two horizontal hues represent two-qubit gates, which 
perform the unitary operation shown on the target qubit if 
and only if the control qubit (indicated by a black dot) is |1) 
(an open circle represents a NOT-gate). Time goes from left 
to right. 
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FIG. 2. Experimental deviation density matrices p^xp for 
l^o) ~ |1)|0)|1), shown in magnitude with the sign of the 
real part (all imaginary components were small), after 2 (a) 
and 28 (b) Grover iterations. The diagonal elements give 
the population difference with respect to the average. The 
off-diagonal elements represent coherences between the basis 
states. Inset: The corresponding ^"'C spectra (^"^C was the 
least signiflcant qubit). The receiver phase and read-out pulse 
are set such that the spectrum be absorptive and positive for 
a spin in |0). 
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FIG. 3. (a). Experimental (error bars) and ideal (circles) 
amplitude of dx^ , with fits (dotted) to guide the eye. Dashed 
line: the signal decay for ^^C due to intrinsic phase random- 
ization or decoherence (for ^^C, T2 « 0.65 s). Solid line: the 
signal strength retained after applying a continuous RF pulse 
of the same cumulative duration per Grover iteration as the 
pulses in the Grover sequence (averaged over the three spins; 
measured up to 4 iterations and then extrapolated) . (b). The 
relative error tr- 



